Effective Lagrangian at Cubic Order in Electromagnetic Fields 
and Vacuum Birefringence 



S. I. Kruglov 

University of Toronto at Scarborough, 
Physical and Environmental Sciences Department, 
1265 Military Trail, Toronto, Ontario, Canada MIC 1A4 

Abstract 

The effective Lagrangian of electromagnetic fields at the cubic or- 
der in field strength is considered. This generalized Lagrangian is mo- 
tivated by electrodynamics on non-commutative spaces. We find the 
canonical and symmetrical energy-momentum tensors and show that 
the vacuum in the model behaves like an anisotropic medium. The 
propagation of a linearly polarized laser beam in the external trans- 
verse magnetic field is investigated. We obtain the dispersion relation 
and refraction indexes for two polarizations. From experimental val- 
ues of the induced ellipticity, reported by PVLAS collaboration, the 
constraint on parameters in the effective Lagrangian is evaluated. 

1 Introduction 

The magnitude of the rotation [1] and ellipticity [2] (reported by the PVLAS 
Collaboration) of a linearly polarized laser beam propagating through a 
transverse magnetic field can not be explained within quantum electrody- 
namics (QED) [3], [1]. This stimulates activities in the theoretical proposals 
on physics beyond the Standard Model. Possible explanations of the phenom- 
ena observed, in the framework of particle physics, are in [5], [6], [TJ, [8], [9], 
|10j . [TT] . [12] . The effect of vacuum birefringence takes place also in Lorentz 
violating electrodynamics [13] . Two popular scenarios, explaining data of the 
PVLAS experiment, include the existence of a new axion-like (spin-0) parti- 
cle (ALP) [6] and/or minicharged particles (MCPs) [8] (see [H], [15] for last 
reviews). The best fit is obtained for MCPs of spin-1/2 [H]. In addition, the 
parameters of ALP are different as compared with the parameters of a QCD 
axion. The case when the coupling and the mass of an axion-like particle 
depend on the temperature and matter density was considered in [16]. This 
can adjust astrophysical bounds to allow for the PVLAS signal. 
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In this letter, we study the induced ellipticity of the laser beam (vac- 
uum birefringence) using the effective Lagrangian at cubic order in the elec- 
tromagnetic field strength which is motivated by electrodynamics on non- 
commutative (NC) spaces. 

Models on NC spaces attract a great interest because NC coordinates 
appear in the superstring theory with the presence of the external background 
magnetic field [17] . In the NC field theories the Lorentz invariance is broken 
due to the fact that the constant parameters 6^ are coupled to tensors, but 
a twisted form of the Lorentz invariance is valid [18]. The NC parameter 
9 is extremely small and astro-physical bounds on Q~ x l 2 are of the order of 
the Planck scale. Experimental bounds on the NC parameters are discussed 
in [11]. Coefficients for the Lorentz violation in electrodynamics within the 
Standard Model Extension (SME) are estimated from Cosmic Microwave 
Background (CMB) radiation in [20]. The Lorentz violation constraints in 
the photon sector from measurements of the linear polarization in gamma-ray 
bursts were considered in [21j. 

We use the Heaviside-Lorentz system of units, and h—c—1. 



2 Lagrangian and Field equations 

It was shown that there is no polarization rotation and ellipticity observed 
in the NC version of electrodynamics [22], [23], [21], [25], [25], W\- The 
effective Lagrangian of the NC version of electrodynamics, at first order in 
the NC parameter 8, is cubic in the electromagnetic field strength. There- 
fore, to explain the data of the PVLAS experiment, we consider the general- 
ized effective Lagrangian of electromagnetic fields at cubic order in the field 
strength. Possible structures, including second-rank "tensor" parameters in 
the Lagrangian, are as follows: 

flW TP TP? fl( 2 ) TP TP TP TP TP? TP TP TP 

u a/3 r a/3 r a/3 T m 1 U T W '> a/3 r a/3 r ^ v , u a /3 r a(3 r fiv r \xvi 

Tp Tp fp ()(") fp TP TP f)(?) fp fp Tp a(&) TP TP TP 

" a /3 r af3 " \u> " fiu , ™ a (i lie*- vf3 r iiv i ™ a (3 l^a v/3^ fivi ™ot/3 V P H u i 

where F^ u = d fl A u — d u A^ is the strength of the electromagnetic field, 
F/J.V = (l/2)^ a /3-F Q/ 3 (^1234 = —i) is a dual tensor. It is easy to verify that 
these structures are not independent and may be converted into only the first 
two terms. So, we have two independent antisymmetric "tensor" -parameters, 
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or four three- "vectors" . As parameters 9^ are not transformed as real ten- 
sors, the Lorentz invariance is broken. Here we investigate the "minimal" 
extension of NC electrodynamics, and therefore, the possible independent 
structure 

(including forth-rank "tensor" -parameters Q^p^v) is not considered here. 

The generalized effective Lagrangian of electromagnetic fields at cubic 
order in the field strength is given by 
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Two constant "tensors" O^p and 9^p are independent. At 9^ = O^p, one 
comes to NC electrodynamics [28J (using Seiberg-Witten map [17]) with the 
accuracy of 0(6 2 ). The Lagrangian (1) can also be rewritten as 

£ = i(E 2 -B 2 )[l + (a-B)-(^.E)]-(E.B)[(/3.E) + ( 7 .B)], (2) 

where the electric field is E{ = iFa and the magnetic induction field being 

Bi = tij k Fj k , a>i = 29f^ — O^, (3i = 6f \ = (l/2)e i j k 6 ( ^ 2 \ and % = i9$ , 



£i = 2i6\ 4 — i6 iA . The parameters Ofy 2 " 1 have the dimension of (length) 2 . 
It follows from Eq.(2) that terms containing parameters c^, (3i violate CP 
- symmetry. The Lagrange-Euler equations lead to equations of motion as 
follows: 



-^9f, {FapFfta) + Oppdp (F a pF ua ) - 9^dp (F^ a F vf3 ) - 0. 
Eq.(3) can be represented as 



(3) 



^-D - V x H = 0, V • D = 0, (4) 

ot 

where the displacement (D) and magnetic (H) fields are given by 



D = — = [l + (a.B)-(e-E)]E-[C9.E) + ( 7 -B)]B 



(E-B)/3-i(E 2 -B 2 )e, 



(5) 
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dC 

H = -= = B[l + (a.B)-(£-E)] + [(0.E) + (7-B)]E 

9B 1 , ( 6 ) 



+(E-B) 7 - - (Er — B J a. 

At the case 0^ = ^l 2 " 1 (a = (£ = 7), we arrive at NC electrodynamics. 
The second pair of Maxwell equations d^F^ = reads 

^-B + VxE = 0, V-B = 0. (7) 

at 

Here we study the propagation of a linearly polarized laser beam in the 
external transverse magnetic field for the case when the effective Lagrangian 
of electromagnetic fields is given by Eq.(l). This consideration generalizes 
the results of [22], [23], [21], [22], [26], [2Z] in the case of two independent 
antisymmetric "tensor" -parameters 9^1 and 9^. 



3 Energy-Momentum Tensor 

Now, we find the energy-momentum tensor to clear up the direction of energy 
propagation. With the help of the standard procedure [29], we obtain the 
gauge-invariant canonical energy-momentum tensor of electromagnetic fields 

T — —T? P f 1 _ IflWfc 1 ^ i lfl(l)p p2 

p p p / p p _i_p p i a( 2 ) p r r ^ ^ 

u [if3 r yv r ~fp \ r p,a r vy 1~ r va r fj.'y ) "a/3 7/3 ^pv^- 

At 0^ = 0^, the canonical tensor (8) converts into one for NC electro- 
dynamics, obtained in [30] . Tensor (8) is non-symmetric, but is conserved, 
dpTfu, = 0. From Eq.(8), we find the energy density £, and the Poynting 
vector P: 

€ = T 44 = [1 + (a ■ B)] — (£ • E) E 2 — (E ■ B) (J3 ■ E) , 

(9) 

T m4 = P = E x H. 

so that the four- vector of the energy-momentum is P M = (P,z£), and the 
continuity equation c^P M = is valid. With the help of Eq.(6), the Poynting 
vector can also be written as 

P = [1 + (a ■ B) - (£ • E)] (E x B) 
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+ - (B 2 - E 2 ) (E x a) + (E ■ B)(E x 7). 

It follows from Eq.(lO) that the direction of the Poynting vector (and the 
energy propagation) is different from the direction of the wave vector k or 
(E x B). So, the vacuum in the model considered, behaves like an anisotropic 
medium (see [26] for a particular case 9^1 = 9^1). 

The symmetric energy-momentum tensor can be obtained by varying the 
action, corresponding to the Lagrangian (1), on the metric tensor g^ v [29J. 
After calculations, we arrive at the symmetric energy-momentum tensor: 

= + -Ol^F^Fpp — 8lpF lfl F p pF lp , (11) 

where the conserved tensor T pu is given by Eq.(8). As the action correspond- 
ing to the Lagrangian (1) is not a scalar, the conservation of the symmetrical 
energy-momentum tensor obtained (11) is questionable. From Eq.(8), (11), 
one can obtain non-zero traces of the canonical and symmetrical energy- 
momentum tensors. For classical electrodynamics, O^p = 9^ = 0, and there- 
fore the trace of the canonical energy-momentum tensor vanishes. It should 
be noted that the modified energy-momentum tensor leads to changing the 
curvature of space-time. This may have an influence on the inflation of the 
universe. 



4 Vacuum Birefringence 

Now we consider the plane electromagnetic wave (e, b) propagating in z- 
direction and perpendicular to the external constant and uniform magnetic 
field B = (B, 0, 0). Then E = e, B = b + B. The rotation of the magnetic 
field, in the PVLAS experiment, does not effect the vacuum birefringence 
within QED calculations [5], [I], and therefore, we consider the stationary 
and uniform external magnetic field. After linearizing Eq.(5),(6) around the 
background magnetic induction field B, one obtains the linearized equations: 

di = EijCj + pijbj, hi = (// ■),,!>! + o,j( j (12) 

where 



l + (a-B) 



S i:j - l, I'>j - !,!),. p i:j = Of'j - ^(7 ■ B) - Btfj, 
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Vi/ = 1 + (ot ■ B) Sij + otiBj + otjBi, Oij 
From Eq.(12) and Maxwell equations 

kidi = k{bi = 0, 

k x e = ub, k x h = 
we find the equation for the electric field e: 



as) 

- J B i £ i + <%( 7 -B)+ 7iJ B i . 



(14) 



-cud, 



k 2 (/i _1 ) w + K (^ 1 ) al k t 6 ib - k 2 (/i" 1 ) 6 ib - h {n~ 1 ) u k 



+UJ 2 €ib + LdEijkkjUkb + ^Pij£jmbkn 



(15) 



e b = 0, 



where is the antisymmetric tensor (£123 = 1)- The homogeneous Eq.(15) 
possesses non-trivial solutions when the determinant of the matrix equals 
zero. To simplify the problem, we consider the case 9±p = 9qj = (£ = 
7 = 0). It should be mentioned that NC field theory preserves unitarity only 
for non-zero space-space non-commutativity, 9 0a = [31], [32] • In addition, 
the bounds on the time-space components 9 0a are much weaker (# -1 / 2 > 
(9(10 GeV)) compared to space-space components (6~ 1 / 2 > 0(10 TeV)) 
[j~9] . Evaluating the determinant for this case (£ = 7 = 0), we obtain the 
dispersion relation: 



A 



A + 2n 2 (a-B) -2(J3-B) 



0. 



where 



A = l + (a-B) 



n 



l + 3(a-B N 



(16) 



(17) 



and n = k/u is the index of refraction. There are two solutions to Eq.(16): 



A = 0, n{ = 1 - 2(q ■ B) 
A + 2n 2 (a-B) -2(/3-B) = 0, 



(18) 



209- B). 



In Eq.(18), we use the expansion in small parameters a, /3. The rij_, n| corre- 
spond to the cases when the electric field of the plane wave e is perpendicular 
(e_LB) and parallel (e||B) to the background magnetic induction field B. So, 
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the speed of light is different for two modes. At the case a = /3, we arrive 
at the result [22] , that the speed of light is shifted equally for both polariza- 
tions. Only at the case a 7^ ft, we have the effect of the induced ellipticity or 
birefringence. In the general case £ 7^ 0, 7 7^ 0, parameters £, 7 contribute 
to indexes of refraction n±, nu. If the angle between the polarization vector 
e and the external magnetic induction field B is 9, then the polarization 
vector at z = is e\ z=0 = E (cos9,sm9) exp(-iut). The components of the 
polarization vector at arbitrary z are given by 

e± = E sin 9 expi (k±z — out) , e\\ — E cos9expi (k\\z — oot^ , (19) 

where k± = n±u, k\\ = nvio. We obtain from Eq.(19) [33] 

a = 6, 6=(k ± -k\\\z=({J3-a)-B y juz, sin 2x = (sin 2a) sin 5. (20) 

One finds from Eq.(20) the induced ellipticity (the ratio of minor to major 
axis of the ellipse) 



1 (/3-a)-B ttL 

* = tanx ^ X ^ -5 sin 26* = '- sin26>, (21) 

2 A 

where u = 2tt/\, A is a wave length. We have used here the smallness of the 
5. As a result, after propagating the distance L, initially linearly polarized 
light becomes elliptically polarized. One obtains to first order in the small 
parameter 5: ip ~ 9 (because the angle of the rotation of the ellipse ip is 
given by tan 2ip = (tan 2a) cos 5). There is no rotation of the polarization 
axis of the ellipse. 

With the help of the preliminary results of the PVLAS experiment [2j 



* = (-3.4 ±0.3) x 10~ 12 ^, L = lm, 

pass 

A = 1064 nm, 0=j, ~B = 5.5 T, e~B = 3.25 x 10~ 10 (MeV) 2 , 
one can find from Eq.(21) the constraint for the parameter difference: 



(22) 



{13 - a)jy ~ 10" 9 (MeV)- 2 , (23) 

where the subscript means the projection on the direction of B. The induced 
ellipticity of the PVLAS experiment can be explained within the effective 
Lagrangian (1). Possibly the improvement of PVLAS dada will change the 
estimation (23) to the lower value. 
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5 Conclusion 



We suggest the effective Lagrangian at the cubic order in the electromagnetic 
field strength which contains two "tensors" Ofy 2 \ This is a generalization of 
NC electrodynamics. At the limit 9^J = 9ffi, one arrives at electrodynam- 
ics on NC spaces (with the help of the Seiberg-Witten map). The Lorentz 
covariance is broken because parameters O^ 2 ^ are not transformed as real 
tensors. Lorentz violating structures at the quadratic order in field strength, 
leading to birefringence in a vacuum without a magnetic field, were discussed 
in [13]. The Lorentz violating operators at quadratic order are constrained by 
astrophysical data [SO] . [2T] , and therefore, we do not include these structure 
in the Lagrangian investigated. 

The density of the energy and momentum, and the canonical and symmet- 
ric energy-momentum tensors are found. The canonical energy-momentum 
tensor is conserved, but the symmetric energy-momentum tensor, obtained 
by varying the action on the metric tensor, is non-conserved. The traces of 
the canonical and symmetric energy-momentum tensors do not equal zero, 
i.e., there is a trace anomaly at the tree level. This anomaly is related to 
the violation of the Lorentz invariance. We show that the propagation of 
the electromagnetic wave in the constant magnetic background and the en- 
ergy propagation have different directions, i.e. the vacuum is similar to an 
anisotropic medium. 

It was proven that the model suggested leads to the induced ellipticity 
which, at the case a = /3, 7 = £ = 0, disappears in accordance with the 
previous results [22] , [23] , [21] , [25] , [26] , [27] . We have calculated the induced 
ellipticity through the parameters a, (3. This relation allows us to obtain the 
constraint on parameters introduced to explain the ellipticity observed in 
the PVLAS experiment. For the case £ 7^ 0, 7 7^ 0, the induced ellipticity 
depends on four "vector" -parameters a, (3, £, 7. 

It should be mentioned that the discussed additions to the Lagrangian 
can not explain all of the PVLAS observations because they do not lead to a 
rotation of the polarization (dichroism). We leave the discussion of bounds on 
parameters introduced coming from astrophysics for further investigations. 
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rotation and ellipticity in vacuum and correct Eq.(22) and our estimation 
(23) to the lower value. 
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